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Abstract: We present a space-time description of hadronization of highly virtual 
quarks originating from a deep-inelastic electron scattering (DIS). Important in- 
gredients of our approach are the time- and energy-dependence of the density of 
energy loss for gluon radiation, the Sudakov's suppression of no radiation, and the 
effect of color transparency, which suppresses final state interaction of the produced 
colorless wave packet. The model is in a good agreement with available data on 
leading hadron production off nucleons and nuclei. The optimal energy range for 
study of the hadronization dynamics with nuclear target is found to be a few tens 
of GeV, particularly energies available in the experiment HERMES. 
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A quark originated from a hard process, converts into colorless hadrons owing to confine- 
ment. Lorentz dilation stretches considerably the duration of this process. While hadrons carry 
a little information about the early stage of hadronization, a nuclear target, as a set of scattering 
centres, allows us to look inside the process at very short times after it starts. The quark-gluon 
system produced in a hard collision, interacts while passing through the nucleus. This can 
yield precious information about the structure of this system and the space-time pattern of 
hadronization. 

The modification of the quark fragmentation function in nuclear matter was considered for 
high-pr hadron production in [|I], for deep- inelastic lepton scattering in p, H, ^j, and for 
hadropro duct ion of leading particles on nuclei in . The data are usually presented in the 
form, 

DeffiZh,PT) ...x 

AD{zh,PT) 

where D{x,pt) and Deff{x,pT) are the quark fragmentation function in vacuum and in a 
nucleus, respectively. 

We treat hadronization of a highly virtual quark, perturbatively a gluon bremsstrahlung and 
the deceleration of the quark as a result of radiative energy loss. We assume that subsequent 
hadronization of the radiated gluons, which includes the nonperturbative stage, does not affect 
the energy loss of the quark. 

The radiation of a gluon takes the time 



tr ^ Tjc^il - a) . (2) 

This expression follows from the form of the energy denominator corresponding to a fluctuation 
of a quark of energy u into a quark and a gluon, having transverse momenta and relative 
shares of the initial light-cone momentum a and 1 — a, respectively. If one calculates radiated 
energy taking into account condition (H) one arrives at the density of energy loss per unit of 
length, which turns out to be energy and time independent ||^ like in the string model. 

In the case of inclusive production of leading particles at Zh I, however, energy con- 
servation forbids the radiation of gluons with energy greater than (1 — Zh)]^. Then, the time 
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dependence of the radiative energy loss can be written as 



AEradit) = f dkl ( 



da av 



dug 



e{l-Zh~a) e{t-tr) 
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where dug/dadk^ = e/ak^ represents the distribution of the number of gluons. The factor 
e = 4:asik^)/37i. 

Ahhough soft hadronization is usually described in terms of the string model, we model 



coupling as{k^) = asik^) at k'^ < k^, in the region which is supposed to be dominated by 
nonperturbative effects, the parameter ko ^ 0.7 GeV is chosen to reproduce the density of 
energy loss for radiation of soft gluons {k^ < k^) corresponding to the string tension, dE/dt = 
) ~ 1 GeV/ fm. This value of ko is consistent with the transverse size of a string corresponding 
to the gluon-gluon correlation radius, Rc ^ l/ko 0.3 fm suggested by QCD lattice results. 

After integrating eq. @) in the soft radiation approximation, we get 



Here we have set ti = (1 — Zh)/xBjfniy and t2 = ti Q^/X^, where xbj is the Bjorken variable. 

Eq. shows that for t < ti, the density of energy loss is constant, dE/dt = —eQ'^/2, 
exactly as in the case with no restriction on the radiated energy 0, § • longer time intervals, 
t > ti more energetic gluons can be radiated and the restriction a < 1 — Zh becomes important. 
As a result, the density of energy loss slows down to dE/dt = — ez/(l — Zh)/t which is a new 
result compared to what was known in the string model. At still longer t > t2, no radiation is 
permitted, but obviously a color charge cannot propagate a long time without radiating which 
must be suppressed by a Sudakov's type formfactor. Assuming a Poisson distribution for the 
number of emitted gluons we get the formfactor, F{t) = exp [—ng{t)], where is the number 
of non radiated gluons. 



In order to calculate a time interval for the leading hadron production (or, better, a colorless 
ejectile which does not loose energy anymore), one needs a model of hadronization and of color 



it by radiation as well, choosing the bottom limit 



in small. We fix the QCD running 




\tiJ 




(5) 
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neutralization. In the large A''^ limit, each radiated gluon can be replaced by a qq pair, and 
the whole system can be treated as a system of color dipoles. It is natural to assume that 
the leading (fastest) hadron originates from a qq dipole made of the leading quark and of the 
antiquark coming from the last emitted gluon. This dipole is to be projected into the hadron 
wave function, \E'(/3, It)-, where (3 and 1 — (3 are the relative shares of the light-cone momentum 
carried by the quarks, and It is the relative transverse momentum of the quarks. The result of 
this projection leads to the fragmentation function of the quark into the hadron, which reads 



dtW{t, u, Zh) , 



(6) 



where W{t, u, Zh) is a distribution function of the leading hadrons over the production time t. 

1 
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Here the quark energy Eq{t) = v — AErad{t)- We have chosen a hadronic wave function in 
the light-cone representation which satisfies the Regge end-point behaviour, |\E'/i(/y, oc 
— + /|.r^/6)^^, where r/j is the charge radius of the hadron. 

Fig.l shows function W{t) for several values of Zh and exhibits the approximate (1 — Zh)i'- 
scaling of the mean production time, tp,. = J dt t W{t), which depends weakly on and 
vanishes at Zh ^ 1- 

Our predictions for the fragmentation function D{zh) depicted in Fig. 2 nicely agrees with 
the EMC data M. 



The production of the leading colorless wave packet with the desired (detected) momentum 
completes the process of hadronization. Any subsequent inelastic interaction is forbidden, 
otherwise a new hadronization process begins and the leading hadron energy falls to lower 
values. Such a restriction means a nuclear suppression of the production rate. 

On the other hand, soft interactions of the leading quark during the hadronization in nuclear 
matter cannot stop or absorb the leading quark |jlO|. Although rescatterings of the quark in 
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Figure 1: Distribution of the hadron production time at u = 30 GeV , 
= 6 GeV^ and = 0.70 ^ 0.95 

Figure 2: Comparison of our prediction for D{zh) with data [13] 
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Figure 3: Comparison of our calculations for the v-dependence of nuclear 
suppression (1) integrated over Zh at = 6 GeV"^ with data [13,14] 

Figure 4: Zh-dependence of nuclear suppression. 



the nucleus result in an additional induced soft radiation, ||10| , |rT|, at the same time the quark 
looses much more energy due to the hard gluon radiation following the deep-inelastic scattering 
just as in vacuum. Thus, the induced soft radiation can be treated as a small correction to the 
energy loss and can be neglected, provided is high enough. 

The transverse size of the colorless wave packet produced in a hard reaction can be small, 
therefore the nuclear suppression is weaker due to color transparency. We take into account the 
evolution of the wave packet during its propagation through the nucleus using the path integral 
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Figure 5: Predictions for the dependence of nuclear transparency at 
z/ = 30 and 75 GeV . Data [I4] at u = 75 GeV are integrated over Zh 

Figure 6: Predictions for the Zh dependence of nuclear transparency at 
zy = 30 GeV vs 
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technique developed in [|T^. Figs. 3-5 show quite a good agreement of our parameter- free 
calculations with available data |13|, O on the z/-, z^- and Q^-dependence of nuclear suppression. 



Unfortunately, there is still no data in the region, Zh > 0.8. 

Note that at high energies many of interesting effects go away or are difficult to observe. 
Nuclear suppression integrated over Zh vanishes (see Fig. 3). The region of high Zh, where 
nuclear suppression is expected to be enhanced (see Fig. 4) squeezes at high z/ towards Zh = I, 
where the cross section vanishes. There is almost no Q^-dependence at high u and moderate 
values of (see Fig. 5). 

We present in Figs. 5 and 6 our predictions for the and z^ dependence of nuclear sup- 
pression for the energy range of the HERMES experiment. We expect the onset of nuclear 
effects at moderate values of (Fig. 5) as well as in the region of Zh > 0.8 (Fig. 6), which can 
be tested by the HERMES experiment. 

To conclude, we have developed a phenomenology of electroproduction of leading hadrons 
on nuclei which is based on the perturbative QCD. Our parameter-free model is in a good 
agreement with available data. We stress that the energy range of the HERMES experiment is 
especially sensitive to the underlying dynamics of hadronization. 
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